An overview on the phase diagram of the frustrated two-leg ladder model 
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Using Density-Matrix Renormalization Group, we investigate the general phase diagram of the 
frustrated two-leg ladder with Heisenberg interactions along legs, rungs and diagonals. We confirm 
that all antiferromagnetic gapped states belong to the same universality class as the Haldane phase. 
In a three-dimensional phase-diagram, we determine a continuous surface with singularities in the 
string-order parameter or its first derivative, corresponding to a transition between two Haldane 
phases with different topological order. Some parts of this transition surface are critical with zero 
gap and vanishing string-order parameter. In the complementary parts, the transition is first order 
with finite gap and string order. The boundary of this surface with the ferromagnetic region is a 
critical end line, when the surface is critical and a triple line anywhere else. Part of this boundary 
coincides with the exactly soluble model proposed by D. V. Dmitriev, V. Ya Krivnov and A. A. 
Ovchinnikov [Phys. Rev. B 56, 5985 (1997)]. 
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I. INTRODUCTION 



A/'-legged spin ladders are formed by assembling TV 
spin-i chains one next to the other. Many weakly- 
coupled ladder systems have now been synthetized. 
Among them, the family Srn-iCun+i02n consists in 
weakly coupled \{n + l)-legged ladders which are ob- 
tained from the Cu02 planes of the parent compound 
SrCu02- After the suggestion that mechanisms sim- 
ilar to those occurring in the Cu02 planes, of cuprate 
ceramics may also lead to superconductivitya, an intense 
theoretical and experimental activity has been developed 
on these systems. The doped spin-ladder compound 
Srii-xCaxCu2iOii is superconducting up to 10 K under 
pressure. Some theoretical tools which are specific to one- 
dimensional systems make the theoretical understanding 
of ladders easier than that of the Cu02 planes responsible 
for the high-temperature superconductivity in cuprate 
ceramics. It is believed that the mechanisms leading to 
superconductivity are similar in both systems and the 
understanding of spin ladders should give some insight 
in the physics of the more complex two-dimensional su- 
perconducting systems. 

A first step in the theoretical understanding of these 
systems is the study of undoped ladders in which the 
spin interactions are of Heisenberg type. The exact so- 
lution of the spin-i Heisenberg chain has been known 
for seventy years from the Bethe ansatza. The quantum 
fluctuations prevent long-range antiferromagnetic order, 
there is no gap in the excitation spectrum and spin-spin 
correlations decay in power law as a function of the dis- 
tance. However, the crossover from ladders to the square 
lattice, by increasing the number of coupled chains, is 
far from being smooth. While ladders with an odd num- 
ber of legs exhibit properties similar to those of a single 
chain, i.e. gapless excitations and a power-law fall-off 
of spin-spin correlations, ladders with an even number 



of legs have a finite energy gap to the lowest spin-1 ex- 
citation and exponential decay of spin-spin correlations. 
These remarkable quantum properties rare remininiscent 
of those, first conjectured by HaldaneEl, concerning the 
spin--| Heisenberg chain with, respectively, n odd and n 
even. Therefore, a natural question arises: how are the 
ladder phases related to the phases in the spin chain? 
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FIG. 1. Frustrated-ladder models: (a) The "Zig-Zag 
chain", (b) Frustrated ladder with equal diagonal interac- 
tions, (c) The more general model studied in this paper 
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For the two-leg ladder, this question has been partly 
answered by WhiteQ. He considered a two-leg ladder with 
additional interaction along one diagonal of each plaque- 
tte [Fig. 1 (a)] and showed by varying this interaction 
that there is a continuous path from the Haldanc to the 
ladder phase which does not cross any phase boundary or 
critical points. Some controversial arguments were pub- 
lished recently by WangO. He considered a ladder with 
additional equal interactions along both diagonals of each 
plaquette [Fig. 1 (b)]. Using Densitg-Matrix Renormal- 
ization Group (DMRG) calculationa3 he proved the ex- 
istence of a first-order transition line and suggested the 
occurrence of a phase-transition from the Haldane-phase 
to a gapped singlet phase with different sensitivity with 
respect to boundary conditions. 

It is therefore interesting to study the general "frus- 
trated two-leg ladder model" [Fig. 1 (c)], including two 
different diagonal interactions in the plaquettes, which 
contains both previous works as particular cases and to 
draw a general phase diagram. We thus consider the fol- 
lowing Hamiltonian: 

L 

H — 2^[ Jj|(Si • Si+2 + S2i-i • S2i+i) + 

i=l 

■/xiS2i • §21+1 + Jx2S2i-l ' 821+2] (1) 

(L is the ladder length and the number of spins is 
N = 2L) and propose a general three-dimensional phase 
diagram for this model. There is an obvious symmetry 
of this Hamiltonian by exchanging Jxi and Jx2 and it 
is more convenient to use: 

S = Jxi + Jx2 and 

V = Jxi- Jx2 (2) 

as parameters. We choose the following energy scaling: 
|2J||| + |J+| + |5|-H|I?| =1 (3) 

Due to the symmetry, it is sufficient to consider V > 0. 
We choose as independent variables x — 2Jjj, y — S, 
z = J± and V is given by the previous relation. The 
variables {x,y,z} obey: \x\ + \y\ + |z| < 1 and are thus 
contained in the regular octahedron represented in Fig. 2. 

Using Density-Matrix Renormalization Group 
(DMRG), Conformal Field Theory and some exact an- 
alytical results, we draw transition surfaces separating 
various phases and identify their critical parts. We review 
the most important results obtained earlier for particular 
cases, including the spin-1 and spin-i Heisenberg chain, 
the frustrated spin-1 and spin-i chains with interactions 
between nearest and next nearest neighbors, the dimer- 
ized chain, and the usual ladder, which are all gathered 
in this general phase diagram. We do not put any re- 
striction on the sign of the four exchange parameters 
and provide a new insight in the regions with negative 
(ferromagnetic) interactions. 
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FIG. 2. The regular octahedron inscribed in a sphere of 
radius 1, representing the parameter-space of our Hamilto- 
nian [Eq. (1)] with the energy scaling of Eq. (3). Some 
remarkable points with their {x, y, z} coordinates are listed: 
the spin-i Heisenberg antiferromagnetic (X,Y,L) or ferro- 
magnetic (X',Y',L') chain, products of independent dimers 
(A,Z) or triplet pairs (A',Z'), the usual ladder (AX and 
ZX segments) etc... The "Zig-Zag" chain represented in 
Fig. 1 (a) corresponds to the surface of the irregular octa- 
hedron (AA',XX',ZZ'). In the gray volume, all interactions 
are negative (ferromagnetic). 
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We begin with the study of the vicinity of particu- 
lar points of this octahedron where the low-lying eigen- 
states can be mapped onto those of the Haldane spin-1 
or spin-i chain (section II), and introduce three differ- 
ent phases: two topologically distinct gapped states and 
the ferromagnetic phase. In our three-dimensional phase 
diagram, the points X, Y, L representing the spin-^ anti- 
ferromagnetic Heisenberg chain appear as particular crit- 
ical points lying in zero-gap critical surfaces separating 
the two gapped Haldane phases. In section III, we briefly 
describe the method that we use to determine transi- 
tion surfaces between different phases: a Density-Matrix 
Renormalization group algorithm, with "suitable" open 
boundary conditions. 

Section IV is devoted to particular planes in the octa- 
hedron corresponding the the ladder model represented 
in Fig. 1 (a). This model can also be viewed as a "Zig- 
Zag" Heisenberg chain with nearest and next-nearest 
neighbor interactions and has already been extensively 
studied. We briefly review the most important results 
obtained for this particular case. Using Lanczos diago- 
nalisation and DMRG, we investigate in section V the 
whole phase diagram in the octahedron and draw transi- 
tion surfaces between the three different phases identified 
in section II. We show that no other phase appears. 



II. IDENTIFICATION OF DIFFERENT PHASES 
OF THE PHASE-DIAGRAM 

A. Evidence for two topologically distinct spin-1 
Haldane phases 

The points Z' (all interactions null, except J± = —1) 
and A' (all interactions null except Jx2 = —1/2) rep- 
resent independent pairs of spins 1/2 coupled ferromag- 
netically along rungs and along plaquette diagonals re- 
spectively. The 3^ degenerate ground-state is identical 
to that of L independent spins-1: Si = (821-1-821-1) at Z' 
and Sj = (S2i+82i-3) at A'. In the neighborhood of both 
points, degenerate perturbation theory is relevant as far 
as the other coupling parameters are much smaller than 
the gap between the 3^ degenerate ground-state and the 
first excited state (i.e. << 1). The split 3^ lowest levels 
can be mapped onto the states of a chain of L interacting 
spins-1 with effective Hamiltonian (at first order in the 
coupling parameters): 



The Hamiltonian (4) for (2J|| -|- Jxi + Jx2) > and 
the Hamiltonian (5) for 2 J|| > and Jxi = [Jxi = is 
satisfied in the planes A'XZ, A'XZ', A'X'Z) represent the 
spin-1 Haldane chain. For Jxi = both spin-1 Hamilto- 
nians (4) and (5) are formally identical. However, they 
represent different topological configurations of the orig- 
inal ladder: the spin-1 operators Si represent the projec- 
tion on the triplet state of the sum of a pair of spins- ^ 
along one rung in Hamiltonian (4) and along one diagonal 
of a plaquette in Hamiltonian (5). 

The most important breakthrough for the understand- 
ing of the antiferromagnetic spin-1 Heisenberg chain 
has beea, the discovery by Affleck-Kennedy-Lieb-Tasaki 
(AKLT)lI of an exactly soluble Hamiltonian differing from 
the Heisenberg model only by addition to the Hamilto- 
nian of a biquadratic term: 



Haklt = ^ 



Si ■ Si4 



(Si-Si+i)2 



(6) 



2 J|| + Jxi + Jx2 



^Si-Si 



+1 



(4) 



Its exact ground-state called "ValencCrBond Solid" 
(VBS) is constructed out of valence bondsB, it is nonde- 
generate and breaks no symmetry. It can be represented 
as schematized in Fig. 3 (a). 

• The spin-1 variables are expressed as a sum of two 
spin-i variables 

• Let "00. /3 denote a spin-1 state in terms of the sym- 
metrized spin-i variables: 

^++ - I + + >, 

iJ+-=^-+ = {\ + ~>+\-+>)|^/2 (7) 

^— = I - - > ■ 



• Each spin-i is contracted with one spin-^ at a 
neighbor site into a singlet state: a "valence bond" . 
Contraction of two spin-i variables into a singlet 

state can be written as e°'^\aj3 > /V2, where e"'^ 
represents the antisymmetric tensor: e++ = e = 
and e+~ = -£"+ = 1 



• The contractions are such that every spin-1 is con- 
nected to both its neighbors by a valence bond and 
the VBS state can be written: 

\^VBS 2-^/Vaifte^^"^^a.fee''^"^ • • • 



in the neighborhood of point Z' and 



in the neighborhood of point A' 



The VBS state, used as variational state for the Heisen- 
berg spin-1 chain, gives an upper bound of the energy 
which is only a few percent higher than the best approx- 
imates through DMRG and is thought to capture the es- 
sentials of the physics of the this model. The VBS state 
has some "hidden" topological long-range orderH. 
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(a) 



^■"tl . . . . . ^tl . . i-tl . . i.-tl . . Sri ...SS-l... it-^^ 



(b) 

(0) (-1) (+1) (0) (0) (-1) (0) (+1) 



FIG. 3. Schematic view of a typical component of the VBS 
state described by Eq. 8. (a) The dashed links represent Va- 
lence-Bonds between two Active neighboring spins- i. (b) Se- 
quence of spin-1 values corresponding to (a); the non-zero 
spins are in a Neel order (-1,1,-1,1,....) 

It is easy to see that each term of the sum in the pre- 
vious relation has 5*^ sequences ordered in a way shown 
schematically in Fig. 3 (b). If we discard the sites with 
S'f — 0, the remaining Sf — ±1 sites have Neel order: 
|1, — 1, 1, —1, 1,-1, • • • >. This topological order is char- 
acterized by the non-local "string-order parameter" : 

0^i\m - l\) = (^S^ exp (^tn "f] S^^ S^^^ (9) 

A normalized string-order parameter: 

d^{d) = '^'^"^^ (10) 

< (5^)2 >2 

has been later introduced by WhiteB. Its value for d ^ oo 
is 1 for the VBS state and 0.84 for the Heisenberg spin-1 
chain. 

From the previous discussion about the neighborhoods 
of points Z' and A' in the phase diagram, it is natural to 
introduce two string-order parameters: 

• OjL{d) obtained with = (Sf^ + Sli_^) in Eq. 9, 
with triplet pairs of spins along rungs 

• Oxid) obtained with = (5|, + S^^^^) in Eq. 9, 
with triplet pairs of spins along one diagonal of each 
plaquette 

These two topologically distinct Haldane phases, that we 
now call HD± and HDx, characterized respectively by 
the string-order parameters 0± and Oi^have been in- 
troduced in a recent work by Kim et al.ll2l, with slightly 
different notations. They are schematized in Fig. 4. Us- 
ing field-theoretical methods, Kim et al. have found 
evidences for some transition lines between these two 
Haldane phases around some particular points of the 



phase diagram corresponding to two-decoupled Heisen- 
berg chains (see sections IV and V). 

With Jxi strictly positive, the full Hamiltonian (5) 
represents the "frustrated" antiferromagnetic Heisenberg 
S=l chain, with interactions Ji = (2J|| -I- J±)/4 and 
J2 — Jxi/4 between first and second neighbors. This 
model has been studied by Kolezhuk et al.til They find a 
sharp discontinuity in the string-order parameter Ox at 
J2/J1 Ri 0.7444, suggesting a first-order transition. The 
significance of this first-order transition will appear in 
section V, when placed in the framework of the general 
phase-diagram. 



B. The antiferromagnetic spin-1/2 Heisenberg chain 

In Fig. 2a, the points X: (2J|[ — 1, all other J's null) 
and Y: {Jxi = Jx2 = 1/2, all other J's null) repre- 
sent two decoupled spin-i antiferromagnetic Heisenberg 
chains, and the point L {J± — Jxi — 1/3, all other J's 
null) represents a single spin-^ antiferromagnetic Heisen- 
berg chain. The exact solutiop of this Hamiltonian is 
known from the Bethc Ansatzcl. There is no gap in the 
excitation spectrum and the spin-spin correlations de- 
crease as power law in terms of the distance. In sections 
IV and V, we shall prove that these points lie in critical 
surfaces separating the two Haldane phases with different 
topological order. 



C. The ferromagnetic state 

In the grey volume A'Y'X'Z' of the octahedron 
[Fig. 2 (a)], all interactions are negative. The ground- 
state is ferromagnetic. It has gapless spin-wave excita- 
tions. The ferromagnetic state obviously extends beyond 
this volume. Its limits will be determined in sections IV 
and V. 

(a) 




(b) 

I OCIIM 

FIG. 4. The two topologically distinct Haldane phases 
HDx (a) and HD± (b). The heavy links represent triplet 
pairing and the dashed links schematize Valence Bonds. 
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(a) 




(b) 



FIG. 5. Growing algorithms used in DMRG calculation 
with free boundaries: (a) in the HDx Haldane Phase, where 
triplet pairing occurs along diagonals, (b) in the HD± Hal- 
dane phase, where triplet pairing occurs along rungs. In either 
case, the ground state is non degenerate. 



III. DMRG WITH SUITABLE BOUNDARY 
CONDITIONS 

Our purpose is to determine transition lines or sur- 
faces separating different phases, in particular the two 
topologically distinct Haldane phases HD± and HDx 
defined in section II. We use the Density-Matrix Renor- 
malization Group (DMRG) method with open boundary 
conditions which gives much more accurate results than 
with periodic boundary conditionaj. However, a well- 
known inconvenience for the Haldane spin-1 chain with 
open boundary conditions is the four-fold degeneracy of 
the ground-state. This degeneracy is easily understood 
through the description in Fig. 3 (a) and Eq. (8) of the 
VBS state, which is thought to capture the essential of 
the physics of the Haldane chain. In Fig. 3 (a), with open 
boundary conditions, there are two spin-i at each end of 
the chain which are not contracted through a valence 
bond to any other half-spin. The interaction between 
these extra half-spins fall off exponentially with L and 
in the long-chain limit, there is a four- fold degeneracy of 
the ground-state due to these two free spins. 

In their DMRG treatment of the spin-1 chain. White 
and Huse remove this inconvenient degeneracy by adding 
at each end an extra spin-^ which forms a singlet with 
half of the last spin-1 in the chaint£l. When applying 
DMRG to the frustrated ladder problem, the adaptation 
of this trick is straightforward and already built in the 
ladder topology. In the HDx phase, we use the bound- 
ary conditions shown in Fig. 5 (a) which are generally 
chosen in treating the ladder problemH. In the HD± 
phase, we use other boundary conditions schematized in 



Fig. 5 (b). If the ground state is viewed as a VBS state, 
in either case Figs. 5 (a) and (b) show that there is no 
"free spin" ("free" means not linked to any another spin 
with a "valence bond") at either end. Consequently, the 
ground-state is non degenerate in the long chain limit. In 
either case our DMRG algorithm for growing the chain 
length proceeds by adding one rung to each of the right 
and left Blocks as shown in Fig. 5. What differs is the 
shape of the initial Block: it is a rectangle in case (a) for 
the HDx phase and a trapezoid in case (b) for the HD± 
phase. 

Using the same algorithm for any point of the phase- 
diagram (most previous DMRG studies on ladders have 
proceeded in such a wayO) generally leads to numerical 
difficulties when crossing the phase-boundary. Suppose 
that we start from a point in the HDx phase with the 
suitable algorithm [Fig. 5 (a)] and follow a path to the 
HD± phase, at the transition to the HD± phase, the 
ground-state, which was non-degenerate becomes four- 
fold degenerate and we need to keep a much larger num- 
ber of states to obtain a reasonable approximation of 
the transition point. To determine numerically many 
transition points in a three-dimensional phase diagram, 
we need a "cheap" algorithm for calculating a transition 
point with a reasonable accuracy in a minimum of com- 
puter time. 

Our algorithm uses the fact that at a phase transition, 
there are some singularities in the first or second deriva- 
tives of the energy. There is a discontinuity in the first 
derivative if the transition is first order. Here, if the tran- 
sition is second order the second derivative is generally 
infinite at the critical point. We thus proceed as follows: 

• We choose a path PoPq going from the HDx phase 
(Pq) to the HD± phase (Pq) and parametrized by 
the variable u. Using the algorithm schematized in 
Fig. 5 (a), we calculate at Pq the energy and its 
first and second derivatives. The first derivative of 
the energy, with respect to x — 2J\\, y — S and 
z = J± is obtained by calculating the correlations: 

qi = < 1 + 4Si • Si+2 > /2 
C±= < 1 + 4S2i-i • S2i > /2 
Cxi= < 1 + 4S2i • S2i+i > /2 (11) 

CX2= < 1 + 4S2i-l • S2i+2 > /2 

The derivative of the energy OEq / du along the path 
is obtained from the preceding correlations. The 
second derivative d'^Eo/dv? is obtained numerically 
through the same calculation with an infinitesimal 
variation on u (chosen at the square root of the 
numerical precision). An approximate second or- 
der polynomial expression of the energy near Pq is 
deduced: 

E = Eo + sdEo/ds+{s^/2)d^Eo/ds^ (12) 
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• Using the other algorithm schematized in Fig. 5 (b), 
a polynomial approximation of the energy around 
Pg is obtained in the same way. 

• A first approximation of the critical point is 
given by the intersection of both polynomials. 

• A new starting point Pi is chosen half way between 
Pq and P(f , a new end point P{ is chosen half way 
between P^ and Pq and the process is reiterated. 

An accuracy of 4 to 5 digits is generally obtained for the 
critical point after 8 to 10 iterations. We keep 30 to 60 
states and use chain length L « 50. Some accuracy tests 
we be given in section V. 



IV. THE "ZIG-ZAG CHAIN" 

The ladder model, with only one diagonal interaction, 
represented in Fig. 1 (a) can also be viewed as a spin-^ 
dimerized next-nearest neighbor Heisenberg chain. Using 
the scaling from Eq. (3), the two-dimensional parame- 
ter space for this model corresponds to the eight faces 
of the irregular octahedron XX'AA'ZZ' (Fig. 6). There 
is an obvious symmetry, in this model, obtained by ex- 
changing rung and diagonal interaction. The symmetry 
line XLX'L'X, corresponding to equal rung and diagonal 
interactions splits the parameter surface into two pieces 
with a one to one correspondence. Due to our non-trivial 
energy scaling [Eq. (3)], this mapping does not corre- 
spond however to simple reflections. 

The AXZ face, in which all interactions are anti- 
ferromagnetic has been studied for many years. On 
the symmetry line XL, the point G corresponding to 
Jj_ = Jxi— ,= 2J|| has been considered by Majumdar 
and Ghoshllj thirty years ago. The ground state is dou- 
bly degenerate and corresponds to a product of singlet 
pairs (or diaiers) along rungs or diagonals. Shastry and 
Sutherlandtfl have later proven that the dimer state along 
the rungs (resp. diagonals) is the exact ground state on 
the whole line GZ (resp. GA). Note that dimer states 
have perfect "string oijder" , characterized by the normal- 
ized order parametersQ Ox = 1 on GZ and 0±_ = 1 on 
GA, respectively. Following the path AGZ, the string 
order parameter jumps abruptly from 1 to at the 
Majumdar-Ghosh point G and 0±^ jumps respectively 
from tp-|Ci|iie. The whole symmetry line LX is a transi- 
tion linctJilJ. According to our notations introduced in 
section II, it corresponds to a phase transition between 
the two symmetrical Haldane phases HDx and HD±. 
Both ends X and L correspond to the spin-i Heisenberg 
chain. The vicinities of these poiala^liave been studied 
through field-theoretical methodsEHlEZlill. The transition 
is first order near X and second order near L. On XL 
the nature of the transition changes from first to second 
order at a tricritical point T. From field-theoretical re- 
sults, at T the spectrum develops a gap between the S=0 



ground state and the first excited S=l states and simul- 
taneously the degeneracy between the two lowest S=0 
and S=l excited states is split. Usin^-jLanczos diagonali- 
sation on finite systems, Emery et alEZI have extrapolated 
in the thermodynamic limit the exchange parameter val- 
ues at which the energy difference between the lowest 
S—0 and S—1 excited states vanishes and found T at 
J\\/J± ~ 0.241. Using the same method, we have found 
no evidence for another tricritical point on the symmetry 
line LX' in the AX'Z face. Hence TLX' is a critical line 
probably up to X'. 

In the A'X'Z' face, all interactions are negative and the 
ground state is ferromagnetic. The ferromagnetic phase 
extends further in the neighboring faces: A'X'Z, A'XZ' 
and AX'Z'. In these faces, the line of transition from 
ferromagnetic to antiferromagnetic state corresponds to 
the exactly solvable Hanailtonian studied by Dmitriev, 
Krivnov and Ovchinnikovllj and hereafter referred as the 
DKO model Hamiltonian: 

Hdko = - J2i=i (S2i-i • - j) 

-i'^-^)Etl (S2i • S2i+1 - i) 

+ ^Eti(Si-Si+2-i) (13) 

The ferromagnetic ground-state of this Hamiltonian is 
degenerate, with a singlet state, which can be written 
explicitlyta. With the scaling defined in Eq. (3), we can 
write the line of transition in the following parametric 
form: 

• A'X'Z face 

x= 2J\={v - I) / {I + 2v - v"^) 

y= S = Jx2 = -y/il + z/2) 

z = J^ = -v{v-l)/{l + 2u-v^) (14) 

with < < 1 

• AX'Z' face 

x= 2J\= {v -l)/{l + 2v -v"^) 

y^S^Jxi = -y{y - i)/(i -^2v~ v'^) 

z= J_L = -j//(l + 2iy- (15) 

with < < 1 

• A'XZ' face 

The line XL' corresponding to Jx2 = J± is a sym- 
metry line. On the upper side of this line the fer- 
romagnetic transition is parametrized by 

X = 2 J|| = {v - l)/(i/2 + 2z/ - 1) 
V= S = Jx2 = -yl{y^ + 2u-l) 
z = J^ = + 2v -I) (16) 

with \ < V <2. On the other side, we have: 
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FIG. 6. The surface of the irregular octahedron 

XX'AA'ZZ' represents the two-dimensional phase diagram of 
the "Zig-Zag chain" model [Fig. 1 (a)]. The symmetry line 
XLX'L'X corresponding to equal rung and diagonal interac- 
tions splits this surface into two parts which map each other 
through the exchange of rung and diagonal. Above (resp. be- 
low) this line, the antiferromagnetic gray (resp. white) area 
corresponds to the HDx (resp. HD±) Haldane phase. The 
black area represents the ferromagnetic state. TLX' and T'X 
are two zero-gapped critical lines corresponding to second or- 
der transition between the two Haldane phases. On TX and 
T'E', the transition is first order. G represents the Majum- 
dar-Ghosh point. From Shastry and Sutherland, on the thin 
dotted segments GA and GZ, the ground-state is exact and 
corresponds to a product of independent dimers 



x= 2J\\ = {v-l)/{2v^ -1) 

y = S=Jx2 = -V{V - l)/(2i^2 _ 1) 

z = J_L = -j//(2//2 - 1) (17) 

with \ <v <2 

In the A'XZ' phase, with ferromagnetic rung and di- 
agonal interactions, the symmetry line XL' meets the 
ferromagnetic region at C (1/7,-2/7,-2/7). From field- 
theoretical resultstj, near X this line is critical with zero 
gap and corresponds to a second order phase transition 
between the two symmetric Haldane phases HD± and 
HDx- Using Lanczos diagonalizations on finite systems 
with N up. to 24 particles and the method previously 
describedLD we have found that this critical line ends at a 
tricritical point T' corresponding to 2 J|| ^ 0.42. From T' 
to C, the transition between the two symmetrical gapped 
phases is first order and C appears as a triple point. 

Figure 6 summarizes the results: 

• The X'LTXT'C part of the symmetry line (equal 
rung and diagonal interactions) corresponds to a 
phase transition between two symmetric Haldane 
phases HD± and HDx corresponding through the 
exchange of rung and diagonal. 

• X'LT and XT' (heavy segments in Fig. 6) are crit- 
ical lines with zero gap and second order transition 
while XT and C'T' corresponds to a first order tran- 
sition with degenerate ground state and finite gap. 

• On either side of the critical line, the HD± and 
HDx phase correspond through a simple symme- 
try. Consequently, the topology of their ground- 
states differ, but they have the same thermody- 
namic properties and same critical exponents near 
the critical line. In the usual sense of "universal- 
ity" refering to critical exponents, they belong to 
the same universality classed. 

Our study of the general ladder model [Fig. 1 (c)] 
will now extend this simple diagram in three dimensions. 
Critical lines will be generalized into critical surface etc... 



V. THE GENERAL PHASE-DIAGRAM 

A. Some exact critical lines and surfaces 

In the preceding section, we have seen that in the sur- 
face of the irregular octahedron XX'AA'ZZ' correspond- 
ing to the "Zig-Zag" chain the ferromagnetic boundary is 
exactly known. We shall extend here this exact bound- 
ary to a part of the whole ferromagnetic volume. We 
have put forward the lines GZ and GA where a dimer 
state is the exact ground state. We shall prove that the 
exact dimer state extends over some planar surface and 
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find exactly the boundary between this surface and the 
ferromagnetic phase. 

Let us rewrite the Hamiltonian [Eq. (1)] as a sum over 
plaquettes: 

H = Ef=i H° with 

Hf = 2J||(Si • Si+2 + S2i--i • S2i+i) + 

J_LS2i-l • S2i + 
2JxiS2i • S2i+1 + 2Jx2S2i-l • S2i+2] (18) 

The eigenvalues of Hf for an isolated plaquette are: 

Ai = -3J_i/2 

A2 = J_i/2 + Jxi + Jx2 

A3= Jx/2-2(Jxi + Jx2) (19) 

A4 = - J±/2 

A5 = -\{Jxi + JX2 + [Jl + 5 Jii 

- 6Jxi Jx2 + 5Ji2 - 2Jx(Jxi + JX2)]^) 
Ae = — \{Jxi + Jx2 — [J±^ + ^Jxi 

- 6JxiJx2 + 5J|:2 - 2Ji(Jxi + Jx2)]^) 

The eigenvalue Ai corresponds to a product of rung 
dimers. It is straightforward to check that it is an ex- 
act eigenstate of the total Hamiltonian H if 2 J|| = S = 
J XI + Jx2- The eigenvalue A2 corresponds to the fer- 
romagnetic state, it is always an exact eigenstate of the 
full Hamiltonian H. A lower bound to the ground-state 
energy per unit length Eq/L of the full Hamiltonian H 
is given by the lowest eigenvalue Xmin of Hf . If the low- 
est eigenvalue is A2, the ground state is ferromagnetic. 
The volume where Xmin — A2 provides an inner envelope 
to the ferromagnetic region. After little algebra, we find 
that this inner envelope is defined by the surfaces: 

w = UxiJx2 + {Jxi + Jx2)(2J|| + J±) + 2J|| = 

(20) 

inside the irregular tetrahedron XX'AA'ZZ' and 

w' = {Jxi + JX2){'2J\\ + Ja_) + 2J\yJ^ = (21) 

outside this tetrahedron. 

On another hand, an outer envelope of the ferromag- 
netic region is provided by a simple calculation of spin- 
waves instabilities in the ferromagnetic jphase. For the 
"Zig-Zag chain" , Kolezhuk and MikeskaEil have already 
pointed out that a simple spin-wave instability calcula- 
tion in the ferromagnetic region gives exactly the DKO 
transition line. Introducing two kinds of bosons for repre- 
senting the spin operators on two different legs, we find 
one optical and one acoustic mode. The energy of the 
acoustic mode at fc ^ is: 

E^^ - -wk^ (22) 

where w corresponds to the expression given in Eq. (20) 
The ferromagnetic state is instable with respect to spin- 
waves for u; > 0. The surface w = thus provides an 



outer envelope to the ferromagnetic region. Inside the 
irregular octahedron XX'AA'ZZ', we have thus provided 
an outer and an inner envelope which coincide. The sur- 
face w — [Eq. (20)] is thus an exact boundary of the 
ferromagnetic region. 

It is also worth paying some particular attention to 
the plane ZBZ'B' of the regular octahedron (see Fig. 2) 
which corresponds to 2J|| = {Jxi + Jx2)- As pointed 
out above, in that plane, the dimer-state with indepen- 
dent dimers along the rungs is an exact eigenstate of the 
Hamiltonian. 

As previously emphasised this dimer state has perfect 
string order Ox = 1 and corresponds to the HDx Hal- 
dane phase. Its exact first order transition line with the 
ferromagnetic phase is given by the equation: 

Ai = A2 i.e. 2J^ + Jxi + Jx2 = (23) 

from point D=(-2/5,-2/5,l/5) to point C = 
(—2/7,-2/7,1/7) corresponding to DKO model (see 
Fig. 7). The other part CZ' of the transition line from 
the ferromagnetic to the HD± phase is given exactly by 
Eq. (20). 

The line ZBZ' with Jxi — Jx2 — J\\ ^0 corresponds 
to the so called "composite-spin model"c3. The Hamil- 
tonian can be written as: 

H = H1+H2 with : 
Hi = 2J|| X^i^Li ' Si+i 
H2 = 2 Jx Ef=i S2i-i • S2i (24) 

where Si = 821-^1 + 821 represents the sum of a pair of 
spins-1/2 on a rung. Hamiltonians Hi and H2 commute. 
The ground state E2 of H2 corresponds to a product of 
dimers along rungs i.e : E°/L = -3J_i/2. It is easy to 
check that it is also an eigenstate of H with eigenvalue 
i?2- The states with all pairs Si = S2i-i + S2i along the 
rungs in a triplet state are degenerate eigenstates of H2 
with eigenvalues E2 = LJ±/2. The subspace spanned 
by these 3^ states corresponds to the spin-1 chain. The 
low-energy spectrum of Hi corresponds to that of the 
spin-1 Haldane chain, and the eigenstates of Hi are also 
eigenstates of H. The ground state of the spin-1 Hal- 
dane chain has an energy Ei « —1.401484039 x 2J|| from 
DMRG calculationsEj and is an eigenstate of H with en- 
ergy El — E^ + E2. Starting from the dimer phase at 
point Z, we have a first order transition to the Haldane 
phase HD± at the crossing between these two energy 
levels El + E2 — E2 corresponding to a critical value 
J±/J\\ « 1.401484039 (point H). The Dimer phase has 
perfect string order Ox = 1 corresponding to the HDx 
phase. At the transition, the order parameter Ox jumps 
from 1 to 0, while 0± jumps from to 0.84, the value 
corresponding to the spin-1 antiferromagnetic chain. 
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this lower bound corresponds to 



Z 




FIG. 7. The (Jxi + Jx2) = 2J|| plane. The dashed Unes 
correspond the intersection of this plane with the "Zig-Zag" 
model and G represents the Majumdar-Ghosh point. The 
black area represents the ferromagnetic region, whose bound- 
aries are exactly known. The line Z'BZ represents the "com- 
posite-spin" model: Jxi ~ Jx2 = J\\- At H, there is a first 
order transition from the Haldane HD± phase to an exact 
Dimer phase which has the same topological long-range order 
as the HDx phase. The heavy dotted line CGH represents 
the first order transition between the Dimer Haldane phase 
HDx (grey area) and the Haldane HD±_ phase. The thin 
dotted line is an exact inner bound to the exact dimer phase. 
C = (-2/7, -2/7, 1/7) is an exact Triple point. 

We have calculated the whole transition line CGH be- 
tween the HD± phase and HDx dimer phase through 
the DMRG scheme explained in section III. The tran- 
sition is first order all along this line and the point C 
= (-2/7,-2/7,1/7) appears as an exact triple point. Note 
that the Majumdar-Ghosh point G belongs to this line 
and is also an exact transition point. From C to G a 
trivial but remarkable exact inner bound to the Dimer 
state is obtained when the dimer energy Ai ceased to be 
the lowest eigenvalue for an isolated plaquette, i.e when 
Ai = Ag (see Eq. 20). With our energy scahng [Eq. (3)], 



J_i = (-2 + 5S + v/8-365 + 4152)/2 (25) 
for 5 > and 



J_i = (-2 - 35 + V8 + 28S + 25S^)/2 (26) 

for 5 < (see the thin dotted line in Fig. 7). Beyond G, 
The upper bound to the dimer phase corresponds to the 
intersection of Ai and A3, i.e. to the segment J± = S. 



B. The general phase diagram 



We shall consider successively the 8 sectors of the tetra- 
hedron XX'YY'ZZ' with the notation [e^,ej,,ez], where 
= ± represents the sign of the coordinate a 



1. The [+,-!-,+] sector, with Jy, J± and S positive 



a. The boundary face XYZ 

The boundary face XYZ, with equal antiferromagnetic 
diagonal interactions Xxi = Jx2 = Jx [Fig. 1 (b)], has 
been studied by WangH. There is an obvious symmetry 
which exchanges J|| and Jx by twisting every other rung 
by an angle tt around the axis of the ladder. ZB is a sym- 
metry line. It is thus sufRcient to study the half trian- 
gle XBZ. From the recent field-theoretical results of Kim 
et al.t2l, near X, which represents two decoupled spin-i 
chains, there is a first order transition between the two 
Haldane phases HDx and HD± along the line J± = S. 
From the previous paragraph there is a first order tran- 
sition between the same phases along the symmetry axe 
ZB at H. Using the DMRG algorithm described in sec- 
tion III, we have calculated the full transition line from 
X to iJ. It is show in Fig. 8. It is in good agreement 
with earlier results obtained by WangO. However, in con- 
trast to the conclusions of Wang, our results show that 
all along this line we have a first-order transition with 
finite gap from the HDx to the HD± Haldane phase, in 
agreement with theoretical-field resultsEJ. 

As pointed out in section III, in contrast to our algo- 
rithm which uses different boundary conditions on ei- 
ther side of the transition line to avoid the degener- 
acy of the ground state, Wang has used everywhere the 
DMRG algorithm with boundary conditions correspond- 
ing to Fig. 5 (a). 
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FIG. 8. The XYZ face with equal diagonal interactions and 
all J's positive. The symmetry axis ZB corresponds to the 
"composite spin" model with a first-order transition from the 
Haldane HD± phase to the dimer HDx phase at H . The 
dashed lines are exact first order transition lines between the 
HD± and HDx phase in the neighborhoods of X and Y (two 
independent spin-i chains) obtained through field-theoretical 
methods). The dots are our first order transition points ob- 
tained through DMRG. They are compared to earlier results 
of Wang (crosses). The behavior of nearest-neighbor pair cor- 
relations, gap and string-order parameter along the lines ZB 
and ZW corresponding to Jy = Jxi + Jx2 are shown in Fig. 9 



As he observed, below the transition line, the ground 
state is degenerate and it is nondegenerate above this 
line. However, his conclusion suggesting a transition from 
a Haldane phase to a "Dimer" phase with different prop- 
erties concerning the behavior with respect to boundary 
conditions is incorrect. If we use everywhere the DMRG 
algorithm with the other boundary condition correspond- 
ing to Fig. 5 (b) we find a degenerate ground state above 
the transition line and a non-degenerate ground state be- 
low! The correct interpretation is, as previously stressed, 
the transition between two topologically distinct Haldane 
phases. Typical sets of curves showing the discontinu- 
ities in the first-neighbor pair correlations [Eq. (11)], gap 
and string-order parameters along the dotted line ZW in 
Fig. 8 corresponding to J|| — Jxi + Jx2 and along the 
symmetry line ZB are shown in Fig. 9. 

b. The OYZ plane 

The OYZ (J|| = 0) plane corresponds to coupled 
dimerized chains. In this plane (Fig. 10) the point L 
corresponds to a simple antiferromagnetic Heisenberg 
chain and Y represents two decoupled Heisenberg chains. 
Strong arguments have beea-given for the existence of a 
critical line joining L to YE3'cfl and a simple analy-tic ap- 
proximation of this line has also been conjecturedcJ. 
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FIG. 9. First order transition between Haldane phases with 
different topological order, (a) Pair correlation C_l along ZW 
(crosses) and ZB (full line); pair correlation Cxi along ZW 
(circles) and ZB (dashed line), (b) Gap along ZW (crosses) 
and ZB (full line), (c) String order parameter Ox along ZW 
(crosses) and ZB (full line); string order parameter 0± along 
ZW (circles) and ZB (dashed line). 
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FIG. 10. The transition line between the two Haldane 
phases with different topological order in the Jy — plane. L 
represents a simple spin-i antilerromagnetic Heisenberg chain 
and Y two decoupled chains. The heavy line YLTi is critical 
with zero gap. On TiL (dotted line), the transition is first 
order. The dash-dotted line is a simple approximation to the 
critical line. 

With our noiitrivial energy scaling [Eq. (3)], it corre- 
sponds to the real root of the third-order equation: 

Jl-S{l-J±f-S^ + 2S^{1-J±)=0 (27) 



We have determined through DMRG the transition be- 
tween the two Haldane phases in this plane (see Fig. 10). 
From X to L, the transition line follows approximately 
the analytic form previously proposed but strongly dif- 
fer beyond L. The critical line extends beyond L up to 
a point Ti that we have determined through the same 
method as that used for the "Zig-Zag" model (cf. sec- 
tion IV). From Ti to CI, the transition is first order with 
finite gap. 

c. The transition surface 

We have exhibited four transition lines, between the 
two Haldane phases with different topological order, in 
four different planes of the phase diagram: XAZ, XYZ, 
OBZ, OYZ. We expect that all four lines belong to a same 
transition surface. Using the DMRG algorithm described 
in section HI, we have investigated this transition surface 
over a fine grid in the {x,y} plane. As a good compro- 
mise, we kept 30 to 60 states in the DMRG scheme and 
chose a 0.02 x 0.02 {x,y} grid. 



X 



FIG. 11. The transition surface separating the two gapped 
Haldane phases with different topological order in the [+,+,+] 
sector. Each point of the grid has been obtained through 
DMRG, keeping 30 states. The dark part of the surface is 
critical with vanishing gap and string-order parameters. The 
clear part corresponds to a first-order transition with finite 
gap and degenerate ground state. The white line which sepa- 
rates both parts of the surface thus appears as a tricritical line. 
The cuts of the surface through the planes XALZ, OBHZ, 
XYZ, and OYZ have boon represented in Figs. 6, 7, 8 and 10 
respectively. The straight line LTGX, joining L (one simple 
Heisenberg chain) to X (two decoupled Heisenberg chains) lies 
in the surface. The Majumdar- Ghosh point G corresponds to 
the intersection of this line with the OBHZ plane. Near X, 
our DMRG results arc consistent with the exact result ob- 
tained from field theory (see text), i.e. the transition surface 
is tangent to the plane J± = S = Jxi + Jx2- 

The results are shown in Fig. 11. The cuts of this sur- 
face through the planes XALZ, OBHZ, XYZ and OYZ 
have already been represented in Figs. 6, 7, 8 and 10 
respectively. The straight line joining L (simple antifer- 
romagnetic spin-i Heisenberg chain) to X (two decoupled 
Heisenberg chains) is a symmetry line in the AXZ plane 
(cf. section IV) and lies in the transition surface. We 
used it for an accuracy test on our DMRG algorithm. 
For the points of our {x,y} grid which lie on the pro- 
jection of this line in the XOY plane, we have plotted 
in Fig. 12 the difference |zf - zf'""'^] between the 
approximate coordinate ^^^^^ ^j^g g^j-face obtained 
through DMRG and the exact value z^'^"'^* correspond- 
ing to LX. 
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FIG. 12. A test of accuracy of our DMRG scheme, with 
only 30 states kept, along the exact transition line LX. The 
accuracy is best (7 digits) near the Majumdar- Ghosh point G. 
It decreases drastically to 3 to 2 digits close to X (two decou- 
pled Heisenberg chains) and also decreases when moving to L 
in the critical regime as the second-order transition becomes 
softer. 

Except in the neighborhood of X: two decoupled spin-^ 
Heisenberg chains, where a much larger number of states 
« 30^ should be retained in the DMRG scheme to obtain 
precise results, we obtain an accuracy which is better 
than three to four digits on Zc- The accuracy is the best 
(7 digits) near the Majumdar- Gosh point G. It decreases 
drastically to 3 ^ 2 digits close to X (two decoupled 
Heisenberg chains) and is of about 5 digits near L, in the 
critical regime. |_. 

Field-theoretical resultsE£l prove that in both planes 
XYZ and XALZ the transition line, near X, is first or- 
der and correspond to 5 = Jj^ . In agreement with these 
exact results, our numerical data indicate that the tran- 
sition surface is tangent at X to the plane 5 = Jj_ . 

Knowing this transition surface, we now determine 
which part of it is critical, with zero gap. We already 
know two intersecting critical lines LTi and YLT. We 
expect that they lie in a same critical part of the sur- 
face which is limited by a line of tricritical points ( "tri- 
critical line") passing through T and Ti. To determine 
this tricritical line, we choose the same methoda_aaJ;hese 
previously used for the "Zig-Zag chain" modelEIlEfl (see 
section IV). 

It is proved from field theory that, when umklapp in- 
teractions become relevant variables, simultaneously the 
spectrum develops a gap between the 8=0 ground state 
and the first excited S—1 states and also the degeneracy 
of the first S=0 excited state and lowest S=l state is split. 
Hence in the thermodynamical limit, the difference 6 be- 
tween the energies of these two lowest excited states is 
zero at the critical value for exchange parameters where 
the gap vanishes. As a simple illustration, let us consider 
the curved line OLY (Figs. 10 and 11) corresponding to 
the cut of the transition surface through the OYZ plane. 



FIG. 13. Critical value of iS = Jxi + Jx2 for which the 
difference 5 between the two lowest excited states vanishes, 
as a function of 1/N, for finite chains with periodic bound- 
ary conditions. We follow the transition line QLY determined 
through DMRG in the OYZ plane (Fig. 10). A quadratic ex- 
trapolation gives {Jxi + Jx2)'^ ~ 0.152 in the thermodynamic 
limit, fixing the position of Ti on the transition line. 



Moving along this curve, we use Lanczos diagonalisa- 
tion on finite systems of length L — N/2, with periodic 
boundary conditions, to calculate the value Sc{N) of the 
parameter S = Jxi + Jx2 at which S vanishes. 
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FIG. 14. Mean square deviation with respect to a 1/A*''^ 
law of the difference between the two first excited levels of 
a finite ladder with periodic boundary conditions. We follow 
the cut of the transition surface through the plane 2J|| = 
determined through DMRG. 
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The results are represented in Fig. 13 as a function 
of 1/7V. A quadratic extrapolation gives Sc{N — > oo) « 
0.152 in the thermodynamic limit. 

A much better accuracy (six digits) on the position of T 
on LX, for the "Zig=^ag chain" model has been obtained 
recently by Eggertlia. He used the field-theoretical result 
that, on LX at the critical point T, the marginal opera- 
tor vanishes and the energy difference 6{N) between the 
two first excited levels is exactly proportional to 
Although we have not been able to generalize this field- 
theoretical result for every point on the tricritical line, we 
have found numerically that except in the vicinity of the 
boundary plane XYZ, this law works remarkably 

well for a large part of the tricritical line. And we have 
used it to improve our accuracy. 
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FIG. 15. The transition surface between the two Haldane 
phases with different topological order in the [—,+,+] sec- 
tor. The dark part of the surface corresponds to the zero-gap 
critical part. All points of the mesh are DMRG results. Our 
DMRG result are not accurate in the vicinity of ferromag- 
netic (X') or antiferromagnetic (Y) decoupled chains. The 
grey areas are simple extrapolations. The critical part of the 
surface is bounded through the tricritical TiT'i white line. 
On TiT'if2 the transition is first order with two degenerate 
ground-states and finite gap. The cut by the plane ALZX' 
corresponds to the "Zig-Zag chain" model (cf. Fig. 6 in sec- 
tion IV). The straight line LX' lies in the critical surface. 



Figure 14 illustrates around the critical value S = 
0.152 found by the previous method (Fig. 13) the mean 
square deviation of S{N) with respect to a 1/A^^ law, 
when we follow the transition line determined through 
DMRG in the OYZ plane. The results are based on 
Lanczos diagonalizations on finite systems with N=16, 
20 and 24 sites. The critical value thus obtained S = 
0.15201 is improved up to five digits. This corresponds 
to the accuracy on the position of the transition line ob- 
tained through DMRG (see Fig. 12). 

Taking successive cuts of the transition surface through 
different vertical planes at J|| = C"*' or iS = C"**' and us- 
ing the same methods, we have determined the whole 
tricritical line TiTY representing the boundary of the 
critical part of the transition surface (heavy white line in 
Fig. 11). Near Y, this line is close to the boundary of the 
surface of transition in the XYZ plane. We do know how- 
ever tha-t it does not intersect it: exact field-theoretical 
resultalj prove that in the XYZ plane the transition is 
first order near Y. 



2. The [-,-!-,+] sector 

The transition surface and its zero-gap critical part is 
determined though the same methods as used in the pre- 
vious subsection. The results are shown in Fig. 15. The 
dark part of the surface is critical with zero gap. It does 
not intersect the X'YZ boundary face and extends fur- 
ther in the [—,+,—] sector. It is bounded through the 
tricritical line TiT'i. On the left of this line ( TiT'ifi) 
part, the transition is first order with degenerate ground 
state and finite gap. The cut through the ALZX' plane 
corresponds to the "Zig-Zag chain" model (see Fig. 6 in 
section IV). The LX' straight line is a symmetry line in 
that plane and it lies in the critical surface. 

3. The [-,-,+] sector 

This sector is particularly interesting since the two Hal- 
dane phases with different topological order compete here 
with the ferromagnetic phase. 

a. The ferromagnetic volume 

From section IV A, the A'X'Z plane corresponding to 
the "Zig-Zag chain" model {Jxi — 0) splits the ferro- 
magnetic boundary in two parts. On the Jxi > side 
the exact boundary of the ferromagnetic domain is the 
surface determined by Eq. (20), taking into account our 
energy scaling [Eq. (3)]. On the other side {Jxi < 0), 
we have calculated the ferromagnetic boundary with an 
accuracy of three to four digits through DMRG (see Fig. 
16). The straight line segment CD corresponding to the 
intersection with the OB'Z plane (cf. Fig. 7) lies in this 
boundary surface. Near X\ we have two weakly cou- 
pled ferromagnetic chains. At first order in perturbation 
theory, our Hamiltonian can be written as: 
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H^ff=LJ\\ + ( +S)S\- S%/2L (28) 



where 



L L-1 

S-=^S2i and S^ = ^S2i+i 



represent the total spins of the two coupled ferromagnetic 

chain. Wc deduce that the plane ( Jj_ +<S = 0) is tangent 
at X' to the ferromagnetic surface. 
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FIG. 16. The surface limiting the ferromagnetic phase in 
the [—,—,+] sector. The line A'CFX' at the intersection of 
this surface with the A'X'Z (Jxi = 0) plane corresponds 
to the exactly soluble DKO model. The front part of the 
surface corresponding to (Jxi > 0) is also known exactly 
[Eqs. (3, 20)]. The back part of the surface corresponding to 
( Jxi < 0) has been determined through DMRG. The straight 
line CD corresponding to its intersection with the OBZ plane 
(see Fig. 7) is exact. The line ACX' corresponding to the 
intersection of the surface with the A'X'Z plane is singular. 
From A' to F, our DMRG results show a discontinuity in the 
gradient to the surface. From F to X, there are singularities 
in higher derivatives. The Fig. 18 show that this line corre- 
sponds to the intersection of the ferromagnetic surface with 
the transition surface between the two Haldane phases with 
different topological order. F is at the crossing of the A'CX' 
line with the tricritical line limiting the critical part of this 
latter transition surface. 
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FIG. 17. The hrst order transition between the HDx and 
ferromagnetic phase in the X'Y'Z face. The dashed lines arc 
the cuts through the exact tangent planes to the ferromag- 
netic boundary at X' and Y' (two independent ferromagnetic 
chains). The dots are our DMRG results, the boundary of 
the black area represents a simple inner bound to the ferro- 
magnetic region obtained through [Eqs. (3, 21)]. 

Similarly, the whole segment A'Y' of the OY' axis cor- 
responds to two decoupled identical ferromagnetic chains 
with alternating (Jxi, Jx2) ferromagnetic bonds. The 
segment A'Y' lies in the ferromagnetic boundary surface 
and near A'Y', at first order in perturbation theory the 
Hamiltonian can be expressed as: 

iJe// =i5-F(2J|| J_l)S^-S^/2L (29) 

where 

L/2 L-1 

ST = IZ(S4i + S4i-3) and S^ = ^(S4i-l+S4i-2) 

i=l i=0 

and the plane (2J|| + J±) = is tangent to the ferromag- 
netic boundary all along the segment A'Y'. 

Although the inner envelope to this boundary surface 
determined by Eqs. (3, 21) is not exact, it is everywhere 
an excellent approximation (see the comparison in the 
X'Y'Z plane represented in Fig. 17). This is due to the 
fact that three different exact transition lines: A'Y', CD 
and A'CXX' belong to this envelope and that the planes 
previously considered arc also tangent to this envelope. 

Prom our DMRG results, it is clear that the line 
A'CrX', at the intersection of the A'X'Z plane ("Zig-Zag 
chain" model) with the ferromagnetic surface and corre- 
sponding to DKO exact model is singular. On this line, 
there are discontinuities in the gradient to the surface 
from A to F, a point whose position will be given be- 
low. From r to X', there are singularities in higher-order 
derivatives. 
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FIG. 18. Transition surfaces in the [—,—,+] sector. The 
dark-grey surface A'X'Y' is the boundary of the ferromagnetic 
phase and corresponds to a first-order transition between the 
Haldane HDx and the ferromagnetic phase. The other sur- 
face corresponds to the transition between the two Haldane 
phases HDx and HD± with different topological order. The 
dark part of this surface is critical with zero gap and string 
order parameters. On the clear part, the transition is first 
order and the gap is finite. This surface abuts the ferromag- 
netic boundary along a line ATX' which lies in the AX'Z 
plane {Jxi = 0) and corresponds corresponds to the exactly 
soluble DKO model. The white line T'lF limiting the critical 
zero gap surface ends on the DKO line at F. The part AT 
of the DKO line is "Triple Line" while FX' is a "Critical-End 
Line". The two heavy lines starting from A' represent the 
plane Jxi ~ 0. 7444(2 Jy -|- J±), where a first order transition 
between two Haldane phases with diff'erent topological order 
has been found by Kolezhuk et al. in the frustrated spin-1 
chain. This plane is tangent to our transition surface. 

b. Triple line, critical- end line 

We have calculated the transition surface between the 
two Haldane phases with different topological order in 
the same way as in the other sectors (see Fig. 18) 

To our numerical accuracy, we observe that the meet- 
ing of this surface with the ferromagnetic boundary oc- 
curs in the X'A'Z plane and corresponds to the DKO 
exactly-soluble model. The X'T'iF part of this surface 
is critical with vanishing gap and string-order parame- 



ters. On the other part ATT'i the transition is first 
order and the antiferromagnetic ground-state is doubly 
degenerated. 

The tricritical line T'lF limiting the critical zero-gap 
surface ends on the DKO line at a particular point r:(- 
0.35,-0.2518,0.1464). According to the usual terminology 
concerning three dimensional phase diagrams, the part 
AT of the DKO line is a line of triple points (i.e. a 
"Triple Line") while FX' is a line of critical-end points 
(i.e. a "Critical- End Line"). 

Our numerical results concerning the coincidence of 
the meeting of the two surfaces with the DKO line are 
enforced by the following theoretical argument: 

• In the region with Jx^ < all interactions except 
J± are ferromagnetic. Consequently singlet pairing 
("valence bond") can only occur along the rungs, 
while triplet pairing can take place along diago- 
nals. Consequently only the HDx Haldane phase 
can appear. We conclude that the transition sur- 
face between the two Haldane phases cannot cross 
the {Jxi = 0) A'X'Z plane. It only exists in the 
Jxi > side. 

• In the Jxi > side, the ferromagnetic boundary is 
exactly known and corresponds to the simple Equa- 
tion 20. This surface has no singularity in its gra- 
dient or higher order derivatives. The meeting of 
the two surface should correspond to a singular- 
ity on the ferromagnetic surface: discontinuity in 
the gradient for the first-order part of the surface, 
or singularities in higher order derivative for the 
critical part. Consequently the transition surface 
between the two Haldane phases do not touch the 
Jxi > part of the ferromagnetic boundary. 

• From the two previous remarks, we conclude that 
if the transition surface between the two Haldane 
phases meets the ferromagnetic surface, it can only 
be in the A'X'Z plane {Jx, > 0) i.e. at the DKO 
line. 

c. The neighborhood of the frustrated spin-1 Heisenberg 
chain 

In section II A we have shown that the neighborhood of 
A' can be described through the frustrated spin-1 chain 
model (Eq. 5) with effective first-neighbor Ji = 2 Jy -I- J± 
and second-neighbor J2 = Jxi interactions. This model 
has recently been studied by Kolezhuk et alEiJ through 
DMRG. At a critical ratio r^ = J2/J1 ~ 0.7444, there 
is a sharp discontinuity in the order parameter Ox sug- 
gesting a first order transition. This critical ratio corre- 
sponds here to the plane Jxi — rc{2J\\ + J±). We have 
represented this plane in Fig. 18. It is tangent in A' to 
our transition surface. Consequently our diagram also 
contains the results on the frustrated spin-1 Heisenberg 
chain, and the transition observed by Kolezhuk et al. 
corresponds to a first-order transition between the HDx 
and HD± phases. 
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4- The [+,—,+] sector 



The transition surface between the two Haldane phases 
in the last of the four sectors with > is represented 
in Fig. 19. Except in the thin tip close to X we have 
found no evidence for a zero-gap critical part. The poor 
accuracy of our DMRG results near X does not allow 
to be more quantitative. Again, near A', our calculated 
surface is tangent to the plane Jxi — 0.7444(2 J||_i- J±) 
corresponding to the results of Kholezhuk et al.tll (see 
above) 

An overall view of the four previous sectors correspond- 
ing to Jj^ > is shown in Fig. 20. 



5. The [—,—,—] sector with with Jy, J± and S negative 



In this sector the phase transitions are between the 
ferromagnetic and the HD± Haldane states. The cor- 
responding surface is known exactly and corresponds to 
Eqs. (3, 20). It is shown in Fig. 21. 



z 




Y' 



FIG. 19. The first order transition surface between the 
two Haldane phases HDx and HD± in the [-!-,—,+] sec- 
tor. The two heavy lines starting from A' represent the plane 
Jxi ~ 0.7444(2 J|| +J±), where a first order transition between 
two Haldane phases with different topological order has been 
found by Kolezhuk et al. in the frustrated spin-1 chain. This 
plane is tangent to our transition surface 



z 




FIG. 20. An overall view of the transition surface in the 
upper half space. The black surface A'X'Y' represents the 
boundary of the ferromagnetic phase. The other surface rep- 
resents the transition between the two Haldane phases HDx 
and HD±. The dark part is critical with zero gap, the clear 
part corresponds to a first-order transition with finite gap. 
Singularities on the lines SIA, f2X, QY, and SIX' are not phys- 
ical: they are related to our normalization condition (Eq. 3). 




FIG. 21. The exact limit of the ferromagnetic phase in the 
[-,-,-] sector 
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6. The [+,—,—] sector 



The boundary face XY'Z' of this sector corresponds to 
the model represented in Fig. 1 (b) with ferromagnetic 
equal diagonal interactions and antiferromagnetic J|| in- 
teractions. There is an obvious symmetry by exchanging 
the parallel and diagonal interaction which is obtained 
by twisting every other rung by an angle tt around the 
longitudinal axis of the ladder. Hence there is a one 
to one correspondence from each point of this face to 
a point of the YX'Z' face in the [—,+,—] sector. The 
two-dimensional phase-diagram of the XY'Z' face is rep- 
resented in Fig. 22. The transition line XE' between the 
two Haldane phases HDx and HD±^ is critical. Near X 
this critical line correspjonds to Jj^ = 5 in agreement with 
field-theoretical resultal3. The critical line remains close 
to this asymptotic result except in the neighborhood of 
the ferromagnetic region, near E'. 




Z' 



FIG. 22. The two-dimensional phase diagram in the bound- 
ary face XY'Z' corresponding to equal diagonal interactions 
[Fig. 1 (b)] . The diagonal interactions are ferromagnetic while 
J|| is antiferromagnetic. There is a symmetry which exchanges 
parallel and diagonal interaction by twisting every other rung 
by an angle tt around the longitudinal axis. Consequently the 
phase diagram of the YX'Z' face is identical. The transition 
line XE' between the two Haldane phases HDx and HD± 
is critical and remains close to its asymptotical limit at X 
( Jx = S) obtained from field theory, except near its meeting 
with the ferromagnetic boundary at E'. 



The three-dimensional phase-diagram of the — , — ] 
sector is shown in Fig. 23. The dark surface represents 
the ferromagnetic boundary seen from inside. It cuts the 
(Jxi = 0) A'XZ' plane at the A'C'Z' line correspond- 
ing to the exactly soluble DKO model. The part of this 
surface which corresponds to {Jxi > 0) is known ex- 
actly from Eq. (20). The other part is calculated through 
DMRG. The other surface represents the transition be- 
tween the two Haldane phases HDx and HD±. The 
clear part corresponds to a first order transition with 
finite gap and string-order parameter. The other part 
is critical with vanishing gap and string-order param- 
eter. The blank tricritical line bounding this critical 
part has been determined through the methods explained 
above. Its intersection with the XOY' plane corresponds 
to 2 J|| R:! 0.75 and the critical surface has a small exten- 
sion in the previous [+,—,+] sector around X. However 
our DMRG scheme is inadequate to determine the crit- 
ical surface so close to X (two decoupled chains) with 
reasonable accuracy (see Fig. 19). 

The symmetry line XC, in the {Jxi = 0) A'XZ' plane 
lies in the transition surface (see Fig. 6b). C represents 
the intersection of the line XL' joining X (two decou- 
pled antiferromagnetic chains) to L' (one ferromagnetic 
chain) with the DKO line. The surface starts at the 
J_L = plane roughly perpendicular to it and it curves 
sharply around XC to meet the boundary plane XY'Z' 
along XE', nearly tangentially to the S = J± plane. 
It abuts the ferromagnetic boundary along the A'C'E' 
curve, which is represented from inside the ferromagnetic 
region by a heavy line in Fig. 23. 

The part A'C of this curve lies in the A'XZ' plane 
and coincides with the DKO curve. Since it corresponds 
to the meeting of two first-order transition surfaces, the 
curve A'C is a triple line. This is justified by the follow- 
ing arguments: 

• The transition surface between the two Haldane 
phases crosses the {Jxi = 0) A'XZ' plane at the 
symmetry line XC. There is a one to one corre- 
spondence on each side of the XC line in the A'XZ' 
plane which maps the HDx to the HD± phase 
by exchanging rungs and diagonal. Consequently 
there is no crossing of the transition surface at any 
other place in the A'XZ' plane. The transition sur- 
face lies on the same side {Jxi > 0) of the A'XCL' 
triangle. 

• On this side of the A'XZ' plane {Jxi > 0), the fer- 
romagnetic boundary is exactly known [Eq. (20)]. 
It has no singularity, which means that the transi- 
tion surface between the two Haldane phases does 
not connect to it. 

• From the two previous remarks, we conclude that 
the transition surface between the two Haldane 
phases abuts the ferromagnetic boundary along the 
curve AC in the A'XZ' plane {Jxi = 0) 
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7. The [+, +, -] sector 



Y 



Z' 




X 



FIG. 23. The 3D phase diagram in the [+, — , — ] sector. 
The cut through the boundary face XY'Z' is shown in Fig. 

22 and the cut through the plane A'XZ' corresponding to the 
"Zig-Zag chain" model has already been studied in Fig. 6 (b). 
The daxk gray surface corresponds to the boundary of the 
ferromagnetic phase seen from inside. For Jxi > 0, the 
ferromagnetic boundary is given exactly through Eq. (20). 
The other side Jxi < is calculated through DMRG. The 
other surface represents the transition between the two Hal- 
dane phases HDx and HD±. The clear part corresponds to a 
first-order transition with finite gap and string-order param- 
eters. The dark part is critical with vanishing gap and order 
parameter. The tricritical white line represents the bound- 
ary of the critical part of the surface. This surface abuts the 
ferromagnetic boundary along the curve A'C'E', which is rep- 
resented from inside the ferromagnetic boundary as a heavy 
line. The part A'C of this curve lies in the A'XZ' plane 
(Jx = 0) and coincides with the DKO line. It is a triple line. 

The other part C'E' of this curve is clearly out of the 
A'XZ' plane. The meeting of this curve with the tricrit- 
ical line is numerically found very close to E'. We have 
not been able however to prove that it coincides exactly 
with E'. We conclude that at least part of the curve C'E' 
is a triple line. 



No phase transition has been found in this sector. The 
ground state is everywhere the HD± phase. 



8. The [-,-!-,-] sector 



Z' 




X' 



FIG. 24. The 3D phase diagram in the [— ,-|-, — ] sector. 
The boundary face YX'Z' is obtained from the symmetrical 
diagram represented in Fig. 22 by changing X to Y and Y' 
to X'. The dark surface represents the boundary of the fer- 
romagnetic phase. Its intersection with the (Jx2 = 0) AX'Z' 
plane corresponds to the DKO exactly-soluble model. The 
iJx2 > 0) part of this surface is known exactly [Eq. (20)]. 
The other part is calculated through DMRG. A part of the 
critical surface (clear-gray grid) corresponding to a second or- 
der transition between the two Haldane phases phases HDx 
and HD± has been determined through DMRG. It is very 
close to the boundary face X'YZ' but cuts it only on YE. A 
rough extrapolation of this surface (array of thin lines) is also 
shown; it indicates that the intersection of the critical surface 
with the ferromagnetic boundary is close to the X'YZ' face 
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The three-dimensional phase diagram of the last 
[—,+,—] sector is shown in Fig. 24. The ferromagnetic 
boundary (dark surface) is obtained in the same way as 
in other sectors: the part corresponding to ( Jx2 < 0) is 
known exactly from Eq. (20). The complementary part 
is calculated through DMRG. 

The cut by the X'YZ' face is obtained from Fig. 22 
(XY' Z' face) through the symmetry which exchanges par- 
allel and diagonal interactions. The transition surface 
between the two Haldane phases HDx and HDx^ lies 
everywhere close to this boundary plane and cuts it on 
YE. 

Only a part of this surface (clear gray grid in Fig. 24) 
has been obtained, using our DMRG scheme, with rea- 
sonable accuracy. It is critical with vanishing gap and 
string order. The array of thin lines in Fig. 24 is a rough 
extrapolation. It indicates that the intersection of this 
surface with the ferromagnetic boundary lies also close 
to the boundary face X'YZ'. 



VI. CONCLUSION 

We have presented an overall view on the three- 
dimensional phase diagram of the frustrated two-leg lad- 
der model. This model contains, as particular cases, the 
spin-1 Heisenberg chain with first and second nearest- 
neighbor interactions, the spin-i "Zig-Zag chain" with 
dimerization and frustration, the usual ladder and a va- 
riety of previously studied ladder models that we have 
briefly reviewed. We have proven that only three different 
phases occur: the ferromagnetic phase and two gapped 
Haldane phases with different symmetry and different 
long-range topological order. In a three dimensional 
phase diagram, these three phases are separated through 
transition surfaces. Some parts of the transition surface 
separating the two gapped Haldane phases are critical 
with vanishing gap and string-order parameters; the com- 
plementary parts corresponds to a first order transition 
with finite gap and two degenerate S=0 ground-states. In 
the two dimensional plane cuts of the three-dimensional 
phase-diagram which correspond to the "Zig-Zag chain 
model" , there is a one to one mapping of these two Hal- 
dane phases through a simple symmetry which exchanges 
rung and diagonal interactions. Consequently, although 
the topology of these two Haldane phases differ, the crit- 
ical behavior concerning their thermodynamic properties 
near the critical surface are the same and they belong 
to the same class of universality. Part of the bound- 
ary of the transition surface between these two Haldane 
phases with the ferromagnetic phase corresponds to the 
exactly-soluble DKO model. We expect that this result 
will suscitate further theoretical work on this model, con- 
cerning in particular remarkable points like F and C on 
the DKO fine (see Figs. 18 and 23). 

Most calculations were performed on the Cray T3E 
from CEA-Grenoble. Numerical data concerning the 



transition surfaces are available upon requesti 
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